A model of transport nonuniversality in thick-film resistors 



m 

O 

o 



C/3 



C. Grimaldi 

Institut de Production et Robotique, LPM, EPFL, CH-1015 Lausanne, Switzerland 

T. Maeder 

Institut de Production et Robotique, LPM, EPFL, CH-1015 Lausanne, Switzerland 
Sensile Technologies SA, PSE, CH-1015 Lausanne, Switzerland 



P. Ryser 

Institut de Production et Robotique, LPM, EPFL, 



CH-1015 Lausanne, Switzerland 



S. Strassler 

Institut de Production et Robotique, LPM, EPFL, CH-1015 Lausanne, Switzerland 
Sensile Technologies SA, PSE, CH-1015 Lausanne, Switzerland 

We propose a model of transport in thick-film resistors which naturally explains the observed 
nonuniversal values of the conductance exponent t extracted in the vicinity of the percolation tran- 
sition. Essential ingredients of the model are the segregated microstructure typical of thick-film 
resistors and tunneling between the conducting grains. Nonuniversality sets in as consequence of 
wide distribution of interparticle tunneling distances. 
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Thick-film resistors (TFRs) are glass-conductor com- 
posites based on Ru02 (but also Bi2Ru207, Pb2Ru206, 
and Ir02) grains mixed and fired with glass powders.^ 
Besides the widespread use of TFRs in pressure and 
force sensor applications,^ their transport properties are 
of great interest also for basic research. The percolating 
nature of transport in TFRs has been reported since long 
time and now it is well documented. '^"^ As shown in Fig.l 
where we reports a selection of previously published data 
on different TFRs,^~* the conducting phase concentra- 
tion X dependence of the conductance G of TFRs follows 
a percolating-like power-law equation of the form: 



(1) 



where Go is a prefactor, Xc is the critical concentration 
below which G vanishes and t is the transport critical 



exponent. 



The values of Go, Xc, and t which best fit 



the experimental data are reported in the inset of Fig.l. 

According to the standard theory of transport in 
isotropic percolating systems, Go and Xc depend on mi- 
croscopic details such as the microstructure and the mean 
value of the junction resistances connecting two neigh- 
bouring conducting sites, while, unless the microscopic 
resistances have a diverging distribution function (see be- 
low), the critical exponent t is universal, i. e., it depends 
only upon the lattice dimensionality D. For D — 3, ran- 
dom resistor network calculations predict t = ~ 2.0,^^ 
in agreement with various granular metal systems, 
or other disordered compounds. 

As it is clear from Fig.l, TFRs have values of t rang- 
ing from its universal limit t ~ 2.0 (filled squares, Ref. 
5) up to very high values like t ^ 5.0 (filled diamonds, 
Ref. 7) or even higher.'^ Despite of their clear percolating 
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FIG. 1. Measured conductances on difi'erent Ru02 (Refs. 
6, 8) and Bi2Ru207 (Refs. 5, 7) TFRs. Solid lines are fits to 
Eq.(l) with fitting values reported in the inset, dashed line 
denotes a power-law with exponent t — 2. 
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behavior, TFRs do not fulfil therefore the hypothesis of 
universality common to many systems and instead belong 
to a different, quite vast, class of materials which display 
nonuniversal transport behavior, that is a regime where 
the transport exponent t depends on microscopic details 
(microstructure etc.). Typical examples of nonuniver- 
sal systems are carbon-black-polymcr composites, and 
materials constituted by insulating regions embedded in 
a continuous conducting phase. ^^^^^ 

Despite that TFRs have been historically among the 
first materials for which transport nonuniversality has 
been reported,^ the microscopic origin of their univer- 
sality breakdown has not been specifically addressed so 
far. In this letter we show that the cross-over between 
universality and nonuniversality reported in Fig.l can be 
explained within a single model whose basic features are 
the peculiar microstructure of TFRs and the tunneling 
processes between conducting grains. 

Before describing our model for TFRs, let us first re- 
call the mathematical requisites for universality break- 
down in random resistor networks. Consider a regular 
lattice of sites and assign to each neighbouring couple of 
sites a bond which has finite conductance g with proba- 
bility p and zero conductance with probability l—p- The 
resulting conductance distribution function is then: 



p{g)=ph{g) + {l-p)6{g), 



(2) 



where S{g) is the Dirac delta-function and h{g) is the 
distribution function of the finite bond conductances. For 
well behaved distribution functions h{g), conductivity is 
universal and follows Eq.(l) with t = to. Instead, as first 
shown by Kogut and Straley,^^ if h{g) has a power law 
divergence for small g of the form: 



lini h{g) oc g~ 

g->0 



(3) 



and a is larger than a critical value ac, then transport 
is no longer universal. Renormalization group analysis 
predicts in fact that 

._ j to if iy + 1/(1 -a) <to . .^ 

^~ \ z/ + l/(l-a) if u + l/{l-a)>to ' ^ ' 

where v ~ 0.88 is the correlation-length exponent for a 

three dimensional lattice. By using to — 2.0 wc obtain 
therefore ac — 0.107. Equations (3) and (4) have been 
shown to arise from a system of insulating spheres em- 
bedded in a continuous conducting material (swiss-cheese 
model), and from a tunneling-percolation model with 
highly fluctuating tunneling distances. Here we show 
that Eq.(3) [and consequently Eq.(4)] arises naturally 
from a simple representation of TFRs in terms of their 
microstructure and elemental transport processes. 

Let us start by considering the highly non- 
homogeneous microstructure typical of TFRs. These sys- 
tems arc constituted by a mixture of large glassy particles 
(typically with size L of order 1-3 /xm) and small conduct- 
ing grains of size $ varying between ~ 10 nm up to ~ 200 



nm. In this situation, the small metallic grains tend to 
occupy the narrow regions between the much larger insu- 
lating zones leading to a filamentary distribution of the 
conducting phase. A classical model to describe such a 
segregation effect was proposed already in the 1970's by 
Pike.'' This model treats the glassy particles as cubes of 
size L $ whose edges are occupied by chains of adja- 
cent metallic spheres of diameter Such chains define 
channels (bonds) which form a cubic lattice spanning the 
whole sample. Let us assume for the moment that a bond 
has probability p of being occupied by a fixed number 
n -|- 1 of spheres and probability 1 — p of being empty. To 
each couple of adjacent spheres we assign an inter-sphere 
conductance 0"^ (i = 1, • • • , n). A random resistor net- 
work can be therefore defined as in Eq.(2) where hn{g) 
is the distribution function of the total channel conduc- 
tance g of n conductances in series: 



n 



(5) 



The high values of piezoresistance (i.e., the strain sensi- 
tivity of transport) typical of TFRs,^ and the low tem- 
perature dependence of transport strongly indicate that 
the main contribution to the overall resistance stems 
from tunneling processes between neighbouring metallic 
grains. Hence, if the centers of two neighbouring metallic 
spheres are separated by a distance r, then the intergrain 
tunneling conductance a is approximatively of the form: 



a = a{r) = aoe-2('-*)/«, 



(6) 



where ctq is a constant which we set equal to the unity, 
^ oc 1/VV is the tunneling factor and V is the inter- 
grain barrier potential. Let us make the quite general 
assumption that the centers of the spheres are set ran- 
domly along the channel, so that the distances r change 
according to the distribution function P{r) for a set of 
impenetrable spheres arranged randomly in a quasi one 
dimensional channel. By following Ref.^^, P{r) can be 
calculated exactly and its explicit expression is: 



P(r) 



-(r-$)/(a„-$) 



e(r - $) 



(7) 



where a„ = (1 + L/n^)'^/2 is the mean inter-sphere dis- 
tance and is the step function. By combining Eq.(6) 
with Eq.(7) the distribution f{a) of the inter-sphere con- 
ductances is then: 

/■OO 

f{a)= / dr P{r) S[a - a{r)] = (1 - a„)a-«" , (8) 
Jo 

where 



e/2 



(9) 



To obtain the distribution function hn{g) of the occu- 
pied channels, we first note that Eq.(5) implyies that g 



2 



is dominated by the minimum inter-sphere conductance 
(Jmin among the set of n conductances in series. Hence 
the small-g Umit of hn{g) is just the distribution function 

/ of CTmin: 



C?Crmin/(fmin) 



, (10) 



which, from Eq.(8) and by setting g ~ amin, leads to: 



hm hn{g) ~ n(l - an)g~ 



(11) 



The conducting bond distribution function behaves 
therefore as Eq.(3) so that for an > etc — 0.107 trans- 
port universality breaks down and t > to- For L = 1 
lim, $ = 10 nm, and ^ = 1 nm this is achieved already 
for n < 90, i. e., slightly less then the maximum number 
L/$ = 100 of spheres which can be accommodated inside 
a channel. 

Our model of universality breakdown in TFRs can be 
readily generalized to describe more realistic situations. 
For example, the number n of spheres inside the occu- 
pied channels can vary according to a given distribution. 
It is also straightforward to rewrite Eq.(7) in order to 
describe cases in which the diameter $ of the spheres is 
not fixed, or to let the size of the insulating grains to 
change by assigning a distribution function for L. It is 
then possible to have different scenarios and, more im- 
portantly, to obtain a crossover from transport universal- 
ity {t = to) to nonuniversality {t > to) within the same 
framework. This reminds the experimental situation re- 
ported for TFRs and summarized in Fig.l. 

Let us now comment on the capability of other ex- 
isting theories to describe transport universality break- 
down in TFRs. At a first glance, the swiss-checsc model 
of Ref. 19 is a natural candidate since the large values 
of L/$ typical of many TFRs may lead to an effective 
continuous conducting phase filling the voids between the 
large glassy grains. However, there are examples in which 
nonuniversality has been reported for TFRs with L/$ 
only of order ~ 5 — 10,^ a value probably too small to 
be compatible with the swiss-cheese picture. Even more 
problematic are the cases for which t = to — 2.0 has 
been measured for TFRs with L/$ ~ 100 (see for ex- 
ample Fig.l), while the swiss-checsc model would have 
predicted t > to- Regarding instead the model proposed 
by Balberg,^*^ for which transport is dominated by ran- 
dom tunneling processes in a percolating network, it is 
important to point out that it was defined by using a 
phenomenological distribution function for the nearest- 
neighbour particle distances very similar to our Eq.(7). 
Balberg argured that such form of P{r) is a reason- 
able compromise between the distribution function of 
spheres randomly placed in three dimensions, and the 
effect of interactions between the conducting and in- 
sulating phases. Instead we have shown that Eq.(7), 
and consequently the power-law divergence of h{g), is 
a straightforward outcome of the quasi one-dimensional 



geometry of the conducting channels in the segregation 
model of TFRs. Despite that our model has been for- 
mulated specifically for TRFs, nevertheless it could be 
applied also to other segregated disordered compounds 
for which tunneling is the main mechanism of transport 
and nonuniversality has been reported. 

In conclusion, we have proposed a simple tunneling- 
percolation model capable of describing the observed 
transport universality breakdown in TFRs. Essential in- 
gredients of the theory are the segregated structure, mod- 
elled by quasi-one dimensional channels occupied ran- 
domly by the conducting particles, and intergrain tun- 
neling taking place within the channels. 

This work is part of TOPNANO 21 project n.5557.2. 
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